In this paper, a moving mesh technique and a non-standard Finite Difference Method are combined, and amoving mesh non-standard Finite Difference (MMNSFD) Method is developed to solve an initial boundary valueproblem involving a quartic nonlinearity that arises in heat transfer with thermal radiation. In this Method, the movingspatial grid is obtained by a simple geometric adaptive algorithm to preserve stability. Moreover, it uses variable timesteps to protect the positivity condition of the solution. The results of this computational technique are comparedwith the corresponding uniform mesh non-standard Finite Difference scheme. The simulations show that the presentedMethod is efficient and applicable, and approximates the solutions well, while because of producing unreal solution, the corresponding uniform mesh non-standard Finite Difference fails.